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Equilibrium states of a closed semiflexible polymer binding to a cylinder are described. This may 
be either by confinement or by constriction. Closed completely bound states are labeled by two 
integers, the number of oscillations, n, and the number of times it winds the cylinder, p , the latter 
being a topological invariant. We examine the behavior of these states as the length of the loop is 
increased by evaluating the energy, the conserved axial torque and the contact force. The ground 
state for a given p is the state with n = 1; a short loop with p = 1 is an elliptic deformation of 
a parallel circle; as its length increases it elongates along the cylinder axis, with two hairpin ends. 
Excited states with n > 2 and p = 1 possess n-fold axial symmetry. Short (long) loops possess 
energies « pEo (nEo), with Eo the energy of a circular loop with same radius as the cylinder; in 
long loops the axial torque vanishes. Confined bound excited states are initially unstable; however, 
above a critical length each n-fold state becomes stable: the folded hairpin cannot be unfolded. The 
ground state for each p is also initially unstable with respect to deformations rotating the loop off 
the surface into the interior. A closed planar elastic curve aligned along the cylinder axis making 
contact with the cylinder on its two sides is identified as the ground state of a confined loop. Exterior 
bound states behave very differently, if free to unbind, as signaled by the reversal in the sign of the 
contact force. If p = 1, all such states are unstable. If p > 2, however, a topological obstruction to 
complete unbinding exists. If the loop is short, the bound state with p = 2 and n = 1 provides a 
stable constriction of the cylinder, partially unbinding as the length is increased. This motif could 
be relevant to an understanding of the process of membrane fission mediated by dynamin rings. 

PACS numbers: 46.70.Hg, 68.47.Pe, 87.10.Pq 


I. INTRODUCTION 

The confinement of semi-flexible polymers or filaments by curved surfaces is increasingly recognized to be a key 
element in a number of physical processes. The best studied examples are presented in biology: DNA, actin filaments, 
and microtubules, for instance, have persistence lengths of 50 nm [I], ~16/im [2] and several millimeters [3j, respec¬ 
tively. Because such polymers need to navigate an intracellular environment crowded by membranes and proteins, 
over and over again one needs to understand how they interact with these structures. DNA in eukaryotic cells, for 
instance, makes a virtue of necessity, wrapping around the core of cylindrical histones to facilitate their condensation 
BH§]; intermediate filaments (such as spectrins or lamins) may adsorb onto membranes and, in the process, modify 
their elastic properties EHS; membrane fission is also frequently driven by the polymerization of tightly winding stiff 
filaments of the protein dynamin, whose winding around the neck of a nascent vesicle is believed to cut the vesicle from 
its parent membrane HMD- Since in these and other examples the persistence length can be quite a bit larger than 
the radii of curvature to which these polymers are confined, it is a good approximation to ignore additional thermal 
fluctuations and focus on ground state solutions, which is what we do in this paper. The opposite limit, confining 
highly flexible filaments, is a classical topic in polymer science, which is reviewed in Ref. m ■ While stiff filaments 
are especially frequent and well studied in the context of biology, they also occur in other situations, and confinement 
issues arise there, too. For instance, single walled carbon nanotubes have a diameter dependent persistence length 
in the tens of micrometer range m , and their adsorption onto surfaces or their interplay with other mesoscopic 
microstructures is currently a subject of considerable interest for engineers as well as physicists mi ns]. 
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The defining feature of semi-flexibility is that the polymer’s persistence length substantially exceeds molecular scales, 
such as the monomer size or thickness. This is physically important, because it means that the energetics of bending 
decouples from the minutiae of the chemical structure and can thus be described to very good accuracy by a continuum- 
elastic Hamiltonian depending largely on geometry. The simplest such energy functional describes a semi-flexible 
polymer as a one-dimensional smooth space curve and quadratically penalizes its Frenet curvature, n (s) mm, 

h b = j J ds k ( s ) 2 , ( 1 ) 

where s is the arc length along the curve and A is the bending rigidity. Since there are no other energy terms 
competing with the bending, from now on we set the bending rigidity to unity, A = 1. A functional variation 5H# = 0 
leads to the Euler-Lagrange (EL) equations which characterize the solutions that minimize this energy. 

Substantially more complicated functionals are not only conceivable but have also been studied, adding, for instance, 
stretching, twisting [2TH25], spontaneous curvature {23] or the higher derivative Frenet torsion |25| . An environmental 
bias may imply a separate dependence on the geodesic and normal curvatures, or even the geodesic torsion [26] . For 
the purpose of this paper we restrict our attention to the elementary energy Eq. Q. The problem is still non-trivial 
because the energy needs to be amended with the constraint confining the space-curve to the surface. One needs to 
accommodate this constraint in the calculus of variations, as described, for example, in m- We only consider the 
limiting case in which the surface is much stiffer than the polymer itself, so that the shape of the surface remains 
unaffected. The opposite limit, in which the polymer is infinitely rigid and the surface adjusts, has recently been 
studied by Bozic et al. |28j , who considered a circular ring constricting an axisymmetric neck. Of course, in many 
physically relevant situations both the polymer and the surface respond to each other’s presence. However, the 
discussion of two elastic objects of different dimensionalities pitting their forces against each other through a mutually 
confining geometric constraint is a long story which will have to await a future treatment. 

The first general discussion of surface-constrained elastic curves a la Eq. 0 was given by Nickerson and Manning 
[291 15(7, who derived the EL equation for this constrained minimization problem. As an application, they studied 
cylindrical confinement with a particular choice of boundary conditions m- Later, Marky and Manning considered 
the wrapping of DNA around a histone octamer, taking into account an adhesive interaction between the polymer 
and its substrate m ■ More recently, the problem was treated by van cler Heijden, who also included a twist degree 
of freedom on isotropic [521 and anisotropic [55 rods with forces applied at the boundaries. Later still, van der 
Heijden et al. [34] studied self-contacts of a cylindrically confined elastic rod. In this work, the elastic curve is treated 
as a Kirchhoff rod or more generally within the context of Cosserat theory. We follow the strategy introduced by 
two of the authors in Ref. E2, where the surface constraint was enforced using a local Lagrange multiplier in the 
variational principle, which can be identified as the external source of stresses to which the polymer is subject due to 
the confinement. This provides direct access to the magnitude of the confining forces, as well as their sign, which will 
vary along the polymer. Formally, the multiplier is identified with the normal projection of the derivative of the force 
vector, which permits one to quantify the loss of Euclidean invariance of the constrained system. In this approach it 
is unnecessary to assume any constitutive relations among the bending moments, the energy density, and the strains: 
the balance of forces and torques follow from the minimization of the energy; these two vectors are also identified. 
For simplicity we restrict our attention to closed polymer loops. 

The equilibrium shapes of semiflexible polymer loops, modeled as elastic curves, confined by a sphere were examined 
in some detail in Ref. [27j (see Refs. [35] 136] for a more physical and realistic treatment, involving numerical analysis, 
see also Eznss] for a detailed study of the dynamics of confined semiflexible polymers employing simulations, as 
well as [52, 00] addressing the effect of electrostatics on the adsorption of polymers rings onto surfaces with spherical 
and cylindrical geometry). The analytical treatment in Ref. [2?J was facilitated by identifying the three conserved 
quantities—corresponding to the three components of the torque—associated with the spherical symmetry. The EL 
equations could be integrated completely in terms of elliptic functions. It was shown that the physics of these states 
depends very sensitively on the length of the loop. Indeed, the ground state itself alternated between two states as 
the length was increased. The forces transmitted to the sphere was also shown to depend both non-monotonically and 
discontinuously on this length. As the length becomes increasingly large, all bound states tend increasingly towards 
geodesic behavior. The aim of the present paper is to investigate systematically the binding of a semi-flexible polymer 
loop to a cylinder, not only by its confinement within it but also by the forced constriction of the cylinder by an 
exterior loop. In the absence of additional adhesive forces the latter has no analogue on a sphere. 

A cylinder is distinguished from a sphere both geometrically (curvature anisotropy) and topologically. Specifically, it is 
non-compact. Each of these features will play a role: some times alone, other times together with unexpected outcomes. 
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We begin by showing how the conservation laws associated with the residual Euclidean symmetries conspire not only 
to yield a quadrature, but also to provide the appropriate observables for characterizing the underlying physics. While 
we lose rotational symmetry about the local surface normal, the underlying axial and translational symmetries still 
provide two conservation laws, one for the axial torque and another for the axial force. These first integrals are each 
of third order in derivatives of the embedding functions. However, the two can be combined into a single equation 
eliminating one derivative, which can be cast as a simple quadrature in terms of the tangent angle. This quadrature 
involves three parameters: the two conserved quantities as well as a Lagrange multiplier associated with the fixed 
length of the polymer. 

We are particularly interested in how the equilibrium states, and in particular candidate ground states, depend on the 
length of the loop. For a given length, such loops are characterized uniquely by two integers, the number of times the 
loop wraps around the cylinder, p , and the number of oscillations it executes in the process (its dihedral symmetry of 
order n). These states differ from their spherical counterparts, also labeled by two integers, in a number of ways. If 
its length is small, the loop does not yet explore the non-compact direction; the reduced symmetry however identifies 
the bound ground state as an elliptical deformation of a circular loop (with n = 1, see Fig. 0a)), a state without any 
equilibrium counterpart on a sphere. As the length increases the loop elongates along the cylinder axis, terminating in 
two hairpin ends. At some point, as we show, the loop self-intersects, thereafter scissoring-in accordion fashion-along 
the axial direction. States with higher dihedral symmetries are unstable initially (a state with n = 2 is illustrated 
in Fig. 0b)). As the length is increased, however, these excited states also grow along the non-compact direction, 
forming folded hairpins. Once they develop, it costs energy to undo these folds. The important point is that beyond 
some critical length, the folded state become snagged energetically on the cylinder; they are stabilized. There is no 
analogous process for loops binding to a sphere. Also, unlike a sphere where states morph periodically into multiple 
coverings of a geodesic circle with increasing loop length, the hairpins that develop diverge increasingly from geodesic 
behavior, while at the same time tending towards a simple characteristic structure. Indeed, the asymptotic energy 
assumes the remarkably simple form, E Ptn —> uEq, independent of p , where Eq is the energy of a circular loop with 
a radius equal to that of the cylinder. It is independent of the topology, so that the energy of these limit states is 
highly degenerate. We describe in detail how the axial torque, the energy and, perhaps most importantly, the forces 
responsible for confinement depend on the length of the loop for each pair of integers, n andp. The direction of the 



(a) (b) (c) 


FIG. 1: (Color online) (a) Elliptic state with n = 1 and p = 1. (b) Two-fold state with n = 2 and p = 1. (c) Folded figure-eight 
state, n = 1 and p = 2. The arrow indicates the local normal vector to the surface, which signifies the direction along which 
confining forces have to act in order to keep the polymer on the surface, with a sign convention further discussed in the text. 

normal force transmitted to the surface obviously depends on which side the polymer is bound to: when the loop is 
inside, positive (negative) values represents outward (inward) forces. If the loop is outside, these sign conventions are 
reversed. The latter are especially interesting because, as we show, an external polymer may constrict the cylinder, a 
happy conspiracy permitted by the topology, an impossibility on a sphere. While in the spherical case the confining 
force on an interior bound loop always pushes it up against the sphere, we see that for cylinders the sign of the confining 
force generally changes along equilibrium bound states. This would indicate that additional forces are required to 
stabilize the state (deriving, for instance, from an adhesion energy) if it is free to unbind. In the case of interior 
bound states, we find that the offending regions are extremely localized in the neighborhood of hairpins, suggesting 
that the state is stabilized if the loop simply takes a shortcut off the surface so as to avoid these regions. If the reader 
feels uncomfortable taking shortcuts, they may assume that an additional constant sticky force is provided where 
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necessary. We also provide an exact description of states confined within the cylinder making partial contact with its 
walls. The corresponding ground state is a planar loop aligned along the cylinder axis making contact on its two sides: 
a vertical safety pin. This state has lower energy than its bound counterpart. In the absence of additional adhesive 
forces, the bound interior ground state-if sufficiently short-is expected to be unstable with respect to deformation 
into the cylinder rotating it into this state. As its length increases, however, it would appear to turn stable with 
respect to such deformations, protected by an ever-increasing energy barrier. A bound exterior state winding the 
cylinder once (p = 1) would, on the other hand, simply peel off relaxing towards a circular loop. If the loop winds 
twice or more times around the cylinder (see Fig. [ljc)), however, this is no longer the case. 1 A figure-eight possesses 
a lower energy than a free doubly-wound circular loop, and would be expected to remain completely attached to the 
cylinder if the loop is sufficiently short. Longer loops will not be completely attached but they would nevertheless be 
expected to embrace the cylinder at least partially, applying a force at the points where they remain in contact. This 
is because the presence of the cylinder obstructs their passage to a singly-wound circular loop. The existence of such 
a motif might be relevant to the constriction of membrane necks, indicating as it does that topology may facilitate 
the constriction process. 


The axial torque also depends in a surprisingly non-trivial way on the length of the loop. It vanishes as this length 
become increasing large in all equilibrium loops. However, it does not do so monotonically: we discover the existence 
of intriguing non-trivial states in which the axial torque vanishes which occur when the length is tuned appropriately. 
These states admit an analytic description. 


The paper is organized as follows: In Sec. |TT] we develop the general framework for axially symmetric binding. The 
first integral of the EL equation, associated with rotational symmetry, is derived in Sec. |II A| This framework is 
applied specifically to the cylinder in Sec. ©■ The equilibrium states of a closed curve are constructed in Sec. |IV[ 


The axial torques, the energies, and the forces they transmit to the surface are determined. This is done perturbatively 
for loop states approximating a multiply-covered circle and numerically otherwise. We end with a discussion and a few 
suggestions for future work in Sec. [V] A number of useful identities and derivations are collected in a set of appendixes 
where the analysis of vertical confined loops, loops with vanishing axial torque and force, and the comparison of elastic 
curves on the cylinder with their planar counterparts are also presented. 


II. CURVES CONSTRAINED TO SURFACES 

Consider a space curve in Euclidean space E 3 , T : s —> Y(s) parametrized by arc-length constrained to lie on a 
surface E. This surface is described in parametric form by the mapping E : (u 1 ,^ 2 ) —> X(u,tt 2 ) £ E 3 . The tangent 
vectors adapted to this parametrization are e a = <9X/<9 u a , a = 1,2, the unit vector normal to the surface is n, the 
induced surface metric is g a t = e a • and the extrinsic curvature tensor is K a b = e a • d&n. The confined curve can 
then also be described as a surface curve T^ : s —X (U 1 (s), U 2 (s)). The surface-bound curve also carries a Darboux 
frame, {T,l,n}, where T = Y' is the unit tangent vector to T, 2 n(C/(s)) is the restriction of the unit normal vector 
along the curve, and 1 = n x T is the conormal tangent to the surface, that is, the vector normal to the curve T but 
tangent to the surface. The structure equations-analogous to the Frenet-Serret (FS) equations-describing how the 
Darboux frame rotates along the curve, are given by [4T 

T' = K g l-K„n, Y = -K g T + T g n, n' = re n T - r g \ , (2) 

where the geodesic and normal curvatures, ti g and as well as the geodesic torsion r g , are defined by 

Kg = T' ■ 1 = l a t b V b t a , K n = —T' • n = t a t b Kat , T g = 1' ■ n = -t a l b I< ab . (3) 

Here t a and l a are the components of the vectors T and 1 with respect to the surface tangent basis e a , a = 1,2: 
T = t a e a , 1 = l a e a . V Q is the covariant derivative compatible with g a t>. Whereas n n and r g depend on the extrinsic 
curvature, ti g is defined intrinsically; it depends only on the surface metric g a b . The FS frame is related to its Darboux 
counterpart through a rotation by an angle w about the tangent direction. This connection between the two frames 
provides a decomposition of the FS curvatures into its intrinsic and extrinsic parts 

tig = K COS LO , tin = K sin UJ , T = T g — <J , (4) 


1 Note that a short interior bound figure-eight would be expected to unwind into the confined planar vertical hairpin state 

2 Here, and elsewhere, a prime represents derivation with respect to the arc-length s. 
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so that K 2 = k 2 + k 2 . The torsion r is given by the difference between the geodesic torsion and u/, the rotation rate 
of one frame with respect to the other. Note that r g involves two derivatives, whereas r involves three. The extra 
derivative is associated with the derivative of ui. 

In equilibrium, one finds that the tension along the curves is balanced by a normal force due to the confinement by 
the surface da 123 , 


F' = -An, 


(5) 


where the tension in the loop F is given by 


F = 


+«n 


C T + ( K 'g + K ^Tg) 1 - (k'„ - KgTg) n . 


( 6 ) 


F represents the tension on a curve segment exerted by the segment with lower value of arclength |25] ; for instance, 
for a spiral traversed upwards, it provides the force exerted on a segment by the region below. The constant C is 
associated with the constraint of fixed length and can be identified with the Hamiltonian [25] , 


The tangential EL derivative £x = F' ■ T vanishes identically, a consequence of the fact that the only relevant degrees 
of freedom are geometrical. 

The projection of Eq. ([5]) onto 1 provides the EL equation 


£i — F / • 1 — n’g + n g 


K 2 g + K 2 n 


~<~ C 


( K l T g)’ 


= 0. 


(7) 


This equation was first derived in Ref. m using a procedure different from that in m- Note that it involves the 
two curvatures as well as the geodesic torsion. In general, the curve will not follow a geodesic with n g = 0. 


The corresponding projection onto n determines the magnitude of the force A transmitted to the surface, 


A = -F' • n = 


-i~ c 


( K s T s) 


( 8 ) 


Thus the normal force is completely determined when the local geometry is known. Its magnitude will vary along the 
contact region even for confinement by a sphere. For a curve confined inside the surface, when A is positive (negative) 
the curve is pushing (pulling) the surface. This expression is missing in the framework presented in [30 . It was 
presented in 


A. Axial symmetry and residual conserved torques 

The integrability of the EL equations for the unconstrained curve is a consequence of the Euclidean invariance of its 
energy, which implies the conservation of both forces and torques. The constrained counterpart will not be integrable, 
for there will not be a sufficient number of conserved quantities; in general, under confinement, one surrenders not 
only translational invariance, but also rotational invariance. The torque about the origin per unit length of the curve, 
M, is given by two contributions: the torque due to the force F about the curve’s position and an intrinsic torque, 

M = Y x F + S , where S = n g n — K n 1. (9) 

Like the tension vector F, the vector M provides the torque on a curve segment exerted by the one with a lower 
value of s. For a free curve the torque vector is conserved, M' = 0, which can be cast in the manifestly translational 
invariant form T x F + S' = 0 [25]. By contrast, for a confined curve one has instead 

M' = ei (Y x 1 ) - A (Y x n) . (10) 

Thus, in a confined equilibrium, M will not generally be conserved, even though £i = 0, because, in general, neither 
the normal force A vanishes, nor are Y and n parallel. 
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Let us now consider a curve confined by a surface with axial symmetry (about the Z axis, say). The axisymmetric 
surface is described by the radial and height coordinates R and Z of the generators, parameterized by the arclength 
£. 3 The azimuthal angle is denoted by ip (see Appendix |a|). The curve is described by the embedding s —>• (£(s), <p(s)), 
or T : s —> Y(s) = R(s)p(s) + Z(s) z, where p(s) = (cos ^(s), sin ip(s), 0). The tangent and conormal vectors are given 
by 


T = cos a ip + sin a e e . 


1 = — sin a ip + cos a e, 


( 11 ) 


where <p(s) = (— sin <p(s), cos <p(s), 0) and e^(s) = R(s) p(s) + Z(s) z; a is the angle that the tangent of the curve 
makes with the azimuthal direction <p. Comparison with the expression (A2| for the embedding functions of an 
axisymmetric surface gives 


sin a 


= t' cos a = Rip'. 


Using Eqs. (Ill and (12), the acceleration T' along the curve is given by 

R r 

T' = (a! — — cot a) 1 — (sin 2 a k± + cos 2 a Ky) n , 

where kj_ and /ty are the principal curvatures along the meridians (generators with constant ip) and the parallels 
(circles of constant £), respectively (see Appendix |A|. Thus the geodesic curvature of T, n g = T' ■ 1, is given by 


( 12 ) 


(13) 


K a = -- 


(R cos a)' 


R sin a 

the corresponding normal curvature n n = — T' • n is given by Euler’s equation for the axisymmetric surface, 


K n = sin 2 a k± + cos 2 a km 


Likewise the geodesic torsion t„ is given by 


T g = sin a cos «(k|| — k_l) . 


(14) 


(15) 


(16) 


Using the above expressions for the Darboux frame of a curve on an axisymmetric surface, along with the expression 
Z = R Ky (see Appendix |a| and the identity cos a Ky = — cos a K n — sin a r g , obtained by a linear combination of Eqs. 


( 151 and ([ 16 |, one finds that the torque about the symmetry axis, M z = M • z, is given by 

rZ = R (— sin a {k +2 K n r g ) + cos a f 9 n — C ) ) + Rn g . 


M z 


(17) 


It is easy to confirm that this component is conserved: projecting Eq. (101 onto z, one gets 

M z ' = M' • z = Y • ( £l l x z - An x z) . 


(18) 


In an axisymmetric surface the three vectors Y, n and z are coplana r an d 1 x z = — sinap — cot aR'ip: therefore, 
M z = — sin aRs\, so it vanishes in equilibrium as claimed. Equation (171 provides a first integral of the EL equation 
|7]), which is not obvious by inspection. 

To reconstruct the curve, the second order differential equation ( fl7| needs to be solved for the angle a with fixed 
M z and C, and appropriate boundary conditions (BCs), which could be specifying the initial angle and its derivative 
(Cauchy BCs) or the angle at the boundaries (Dirichlet BCs). Once a is known, the position on the surface is 
determined using the relations (12). The two constants M z and C can be tuned to fix the length of the curve, L, and 
total azimuthal angle turned, Aip. 


III. CYLINDRICAL CONFINEMENT 

Now we apply this framework to treat curves lying on a cylinder. Despite the apparent simplicity of the cylindrical 
geometry, the determination of the equilibrium configurations of curves on a cylinder is not straightforward. However, 
as explained below, the additional translational symmetry facilitates the integration of the EL equation. 


3 Derivatives with respect to i are denoted by an overdot ,' = d{. 
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A. Derivation of the quadrature 


The height function is given by the meridian arc length, Z = t. We normalize all lengths in terms of the cylinder 
radius Rq. Scaled quantities are denoted by lowercase letters or by an overbar; e.g., the scaled height function is 
z := Z/ Rq. Curves on the cylinder are parameterized by scaled arc length s/Rq . 4 


The scaled geodesic curvature is given by 


Kn := Rq K a = ot . 


(19) 


The two scaled principal curvatures are k± := Rq k± =0 and Ky := R 0 Ky = 1, so that the scaled normal curvature 
and geodesic torsion are given by 


n n := Rq k u = cos 2 a ; r g := Rq t 9 = sin a cos a . 
Inserting these expressions into the EL equation ([T]) , it is easily seen to read (c = Rq C) 


CSC 


n 3 


e\ = a'" + a' ( -Z- + ^ cos 4 a — 6 sin 2 a cos 2 a — c) = 0 . 


The scaled first integral, Eq. (p~7|) (m := Rq M), reads 


to := m • z = — sin a yn' g + 2n n r g ) + 


cos a 




( 20 ) 


( 21 ) 


( 22 ) 


One could now solve this third order differential equation (second order for a) following the procedure outlined at the 
end of the previous section. However, it is possible to do better. We have yet to exploit the translational invariance 
along the cylindrical axis. When we do, we obtain a quadrature for a. 

Translational invariance along z implies that the projection of the scaled stress vector f := Rq F along the cylinder 
axis is constant. Projecting f onto z, and making use of definitions of T and 1 given in Eq. (Ill, one identifies the 
first integral, 


f := f • z = cos a (k' 9 + 2K n r g ) + sin 
It is straightforward to confirm that /' = Rq cosaei. 


a 


(23) 


Both first integrals are of second order in derivatives of a (this dependence enters through n' g in Eqs. (22) and (23)). 
One can now take an appropriate linear combination to eliminate this second derivative: specifically, the definition 


k 2 — K 2 

//(a) := to cos a + f sin a = y -- — c, 

will lead to a quadrature for a. Substituting for n g and n n one gets: 

l cos^ rv 

-(a 1 ) 2 + U(a) = c, where U(a) = -—-/u(a). 


(24) 


(25) 


It is not obvious at the level of the EL equation that two integrations are possible. However, it should be remarked 
that the poor man’s derivation of the quadrature using a variational principle adapted to symmetry demystifies this 
“coincidence”. It is provided in Appendix [B] a special instance of a more general construction presented in [27[. 

One may ask if any additional information is to be gleaned by taking some other linear combination of / and to. 
Consider 


ic(a) := —to sin a + f cos a = n' g + 2 n n r g . 


( 26 ) 


4 Since there is no possible confusion, the scaled arc-length is also denoted by s. 










In terms of a this second order differential equation reads 


a" + 2 cos 3 a sin a — v(a) = 0 . 


(27) 


On the one hand, taking into account expression (191, one has // = v n a . On the other hand, using the identity for 
the cylinder R' n = —2R g f g , the derivative of the righffhand side of Eq. (241 can be written as 


K 2 - K 2 
K g K n 


— C 


= («/ +2(t„r s )K ; 


9 * 


(28) 


These two relations confirm that Eq. (26) follows from differentiation of Eq. (241. Alternatively, substituting identity 


(281 in Eqs. (221 and (231) for m and /, and integrating them, one obtains the quadrature (241. Further differentiation 
of Eq. 

only reproduces Eq. 


(|26|) and substitution of the EL equation to replace K g in favor of terms quadratic in the Darboux curvatures 


(241. 


The quadrature (251 provides us with an analogue of a particle with unit mass and energy c in a potential (7(a). 


While the hard work lies ahead, this will prove to be very useful throughout this paper. That c plays the role of the 
energy in this identification is particularly appropriate, because c is also the conserved quantity, i.e., the Hamiltonian 
function, associated with the absence of an explicit dependence on s in the bending energy density, as detailed in 
Appendix [B] 

The first term in (7, proportional to cos 4 a (or k^), reflects the anisotropy of the sectional curvatures on the cylinder, 
which breaks the rotational symmetry about the normal vector. As a result, despite the fact that the cylinder is 
isometric to a plane, its extrinsic geometry is anisotropic and elastic curves on a cylinder behave differently from their 
counterparts on the plane. It is precisely this term that distinguishes elastic curves on the cylinder from their planar 
Euler elastic counterparts, which are integrable in terms of elliptic functions, as described in Appendix [Cj This can 
be seen by introducing the change of variables a — >• if) + [3, f — > Fsm/3 and m —>• F cos/3. Now one has p —> F cost/) 
and v — ► —F simp, so without the term and its derivative, the quadrature (251 and the first integral (|26[ 
to Eqs. 


(C8). 


reduce 

If the curve lies close to a parallel or the cylinder has a very large radius, R 0 1, so that n n —> 0, 
0 and K g becomes the FS curvature of the curve on the plane, the two descriptions coincide. The comparison 
between the two is discussed in Sec. IIV B1 


The two coordinate functions of the curve on the cylinder, are determined by the relations (121, which in terms of the 
scaled quantities reduce to 


Z = sin a , p' = cos a . 


(29) 


B. Total bending energy and constraining force 


The scaled bending energy of the constrained curve is given in terms of the angle a by 

1 


hs := Ro Hb = 


J ds (a' 2 + cos 4 a). 


Using the quadrature (241, the derivative can be eliminated in favor of a, so that 

liB = J ds (cos 4 a + / sin a + m cos a + c) . 

Taking into account the relations z' = sin a and ip' = cos a, the last terms can be integrated 

hs = / d s cos 4 a + fAz + mAp + cl , 


(30) 


(31) 


(32) 


where Az = AZ/Rq. Ap and l = L/Rq are the total scaled height, winding and scaled length of the curve. For a 
closed curve A z = 0 and A<^ = 2-7T p. This relation makes explicit the respective interpretations of /, m, and c as 
coupling to the height, rotational extent and length, respectively. 
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Likewise, using the quadrature (251 and Eq. (27) to eliminate derivatives of a in the expression ([8]) for the magnitude 
of the scaled local constraining force on the cylinder, A = -RqA, can be expressed as 

A = 2 S 2 (5 — 6 R n ) + (6 — 11 R n ) p, — 4 f 2 + 6 c (1 — 2R n ) , (33) 

or in terms of a and constants m, f and c, 

A = 2 cos 4 a (5 — 6 cos 2 a) — 6 c cos 2a — 3 / sin a (5 cos 2 a — 2) — 5 m cos a (3 cos 2 a — 2) . (34) 

It depends only on the local value of the turning angle. 

IV. CLOSED EQUILIBRIUM CONFIGURATIONS 


Now we focus on isolated loops, i.e., closed curves without external sources of tension other than those due to 
confinement. 5 Let the loop have length L = 2ttR, p be the number of times it winds the cylinder and n the number of 
periods in completing a circuit of the loop. Loops are characterized by p and n: loops with different p are homotopically 
inequivalent; n distinguishes equilibrium states within a given homotopy class. We denote the sequence of loops with 
fixed values of these two integers, generated by varying the length of the loop by L p n . Define the scaled total length 
l = L/Rq = 2tt r, where r = R/Rq is the scaled radius of the loop. Since l > 2np, the scaled radius is bounded from 
below by the number of windings, r > p. The length is expressed in terms of the scaled excess radius A r = r — p > 0. 

We examine loops that oscillate symmetrically about a parallel, which we set a,t z = 0 and refer to as the equator. 
This implies that the odd term linear in sin a vanishes, so / = 0. Note that the potential U appearing in (251 


now possesses up-down symmetry, a —> —a. However, the constant m does not necessarily vanish which breaks the 
left-right symmetry a —> n — a. 


The tangent angle a will be bounded in a closed loop. It will now oscillate symmetrically between the turning points 
of the even potential, and Q-m- These turning points will occur at the intersection of the loop with the equator; 

a = 0 occurs on the two extremal parallels. 


The two parameters m and c are determined by imposing the periodicity on the coordinates, required by the closure 
constraint. Note, first of all, that closure in the axial direction is consistent with / = 0. To see this, note that for a 
closed curve, after one half period one has Ax = 0 (T8J 20]; using the quadrature (251 and the relation ([29] connecting 
z to a one has 


1 


sma 


A z = —= / da __ 

a/2 J-O.M y/ c ~U{a) 


= 0 . 


(35) 


This is because U(a) is an even function of a so that the integrand is odd. 


The remaining two parameters m and c can be related to the three geometric quantities p , n and Ar by integrating 
the quadrature (251 along with the relation connecting to a, (Eq. (291) using the fact that the azimuthal range A ip 


is quantized; i.e., A</?/(27 t) = p. One obtains 


Ar = 


p = 


An f aM 
\[2 2 ir Jo 
An r M 


da 

da 


1 


s/ c — U (a) 

COS O' 


\f22ix Jo yjc - U (a) ' 


(36a) 

(36b) 


In Eq. (36b) we have used the fact that a completes n periods while the loop completes p trips around the cylinder. 


Before analyzing the loops in the nonlinear regime, we first solve the quadrature in a perturbative manner for loops 
which deviate slightly from a circular loop covering the equator p times. The insight we obtain will prove useful when 
solving the quadrature (25) numerically. 


5 Recall that we ignore polymer twist, so one might wonder whether this could change the solutions. However, if the polymer is rotationally 
symmetric, having no “preferred side” for it to adhere onto the surface and no preferred bending direction, any local twist does not 
couple to the surface constraint, and so we only need to consider the global constraint. If we now start with an untwisted circle and 
imagine extending its length in such a way that no twist is added, ignoring twist right from the start is permissible. 
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A. Perturbative solutions about a circular loop 


To examine deformations of a circular loop perturbatively, we let oq represent a small perturbation about ao = 0. 
We set / = 0 and expand the constants, to = Too + rn 2 + ... and c = Co + C2 + ... , 6 . The subindices represent the 
order in the expansion. The approximation of the quadrature (251 at lowest order reads 


1 


to 0 = -c 0 


(37) 


Note that in o and Co appear in the combination too + Cq which is constant independent of n and p. At second order 
the quadrature describes harmonic oscillations, 


where we have defined 


\ (a'i) 2 +V 2 (ai) = E 2 , 


V 2 {a\) = a 0 , o 0 = to 0 /2 + 1 and E 2 = m 2 + c 2 . 


(38) 


(39) 


Thus, the behavior of the perturbation can be described in analogy to a particle of unit mass and energy E 2 moving 
in a quadratic potential V2. The parameter ao is necessarily positive. Otherwise, V 2 would be a downward parabola, 
a would not be periodic, and the corresponding loops would be open. 


The two turning points of V 2 , corresponding to the maximum and minimum values of a, are ±aMi> where 

V^E~2 


OtM 1 = 


and q := 2ao = Too + 2 . 


(40) 


Solving Eq. (38) one identifies 


a 1 = cslm 1 cos q s ; 


combining Eqs. (371 and (401, Too and Cq can be expressed in terms of q 

2 „ 3 2 

m 0 = q — 2 , c 0 = -- q . 


(41) 


(42) 


The condition that the periodic perturbation should complete n periods, while the curve winds p times around the 
cylinder implies that the wave number is given by the ratio of periods to windings, 


q = 


P 


(43) 


Notice that too does not vanish because q is a rational number. A critical axial torque is required to deform a loop 
of radius Rq into a state with a given ratio q. In particular, too is negative when n = 1 and positive for all higher 
values of n. This sign change has physical significance. Superficially, states with the same p and different values of n 
may appear to be identical when Ar = 0. However, as we will see, the different values of Too and cq in these states 
determine the critical buckling forces required to initiate the deformation of the circular loop into the appropriate 
state. 


From the relation (401 between E 2 and q?m 1, one concludes that the sum of their second order corrections is constrained 
by the amplitude 


rp , r 2 

E 2 = m 2 + c 2 = — a M !. 


(44) 


First order terms are not considered because they will be proportional to first order changes in the length of the circle, which vanishes 
on account of the geodesic character of the parallel. 
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To identify a second relation between m 2 and C 2 , we use the conditions (361, which relate the excess radius and the 
azimuthal range in terms of the potential. Since to first order the turning points are !«« 1 , the denominator in these 
expressions can be expanded as 


c — U(a) ss (a 2 M1 - aj){x - Zaf), 


(45) 


where 


X=\{<1 2 + m 2 ) 1 , £=^(l8 + g 2 ) 


thus, 


\/2 \/c— U(a) q\Ja 2 M 1 - af [ q 


1 + ^2 ( a M 1 + a i ) y ) 


Substituting this expression into the two conditions (361 and integrating, one obtains 


Ar=^ (3(a 2 M1 -m 2 ) , m 2 =^-|-^a : 


2 

M 1 ■ 


(46) 


(47) 


(48) 


From these two equations one finds that the excess radius is proportional to the square of the perturbation amplitude, 

(49) 


A P 2 


Using this relation in Eqs. (441 and (481, one can express m 2 and c 2 in terms of the excess radius: 

m 2 = - ( 3 - ]- q 2 ^] Ar, c 2 = - (J-q 2 - 9] Ar. 


P 


P 


(50) 


To complete the construction of the constrained curve, one needs to determine its position coordinates £ and ip by 
integrating the relations (291. To first order, 

A = <* 1 , ip[ = 0 ; (51) 

consequently, this confirms to first order that a non-vanishing value of the conserved axial force / is inconsistent with 


closed deformations of a parallel circle. Note first that, if / ^ 0, the potential V 2 defined in Eq. ([59 1 gains a linear 

osure: 


term so that z\ picks up a helical contribution proportional to arclength, which is inconsistent with c 

Z\ = — \ fp~Kr sin qs , ip 1 = 0 . 
n 


(52) 


In contrast to its spherically confined counterpart, the elliptic deformation of an elastic loop on a cylinder with n = 1 
and p = 1 is not identified with a rotation, but with a tilt of the loop extending its length. 


1. Energy and confining force of perturbed loops 

In the quadratic approximation the bending energy hs is given by 

hs ~ tt (p + (2 q 2 - 3) Ar) . (53) 

Using the fact that r _1 « p — Ar, the energy can also be expressed in terms of the energy of the circle with reduced 
radius r, h c = ir/r, as 


h B « he [1 + 2 (q 2 - 1) Ar] . (54) 

It is now manifest that at lowest order the energy of the deformed loop is given by the energy of the circular loop, 
h B ~ h c . The second order correction to the energy of the loop vanishes for elliptic perturbations with q = 1, whereas 
it increases (decreases) if q > 1 (q < 1). 
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To second order, the magnitude of the normal force A given by Eq. (341 reads 

1 


A » q z — 1 + 


P 


q (11 cos 2 qip — - ) + 7 — 2 cos 2 qip 


A r. 


(55) 


For the elliptical confined states with q = 1, only the second order correction is non-zero, A « 18 (cos 2 ip — |) A r. It 
vanishes in the limit A r —> 0. Note that A is not positive everywhere, assuming negative values in the two intervals 
ip £ (7t/3, 2/37t) and ip £ (— 2/37T, —7t/3), describing the neighborhood of the turning points where the height is a 
maximum and a minimum. On account of this, if the loop is bound inside the cylinder and free to detach, one expects 
that it will tend to rotate into the vertical plane, touching the cylinder tangentially at two points. Such states are 
described exactly by two segments of a planar Euler elastic curve as described in Appendix [Cj A loop bound to the 
outer side of the cylinder, free to detach, will tend to detach on the equator first before reforming as a circular loop. 


Note that a constant normal force Ao = q 2 — 1 persists in the limit Ar —► 0 in all other states. This behavior 
is analogous to the Euler instability associated with n-folds confined by a sphere m 7 8 In general, Ao is positive 
(negative) if q > 1 (q < 1), so initially the loop will push (constrict) everywhere on the cylinder. In particular, 
an interior bound curve with p = 1 and n-fold dihedral symmetry will push on the cylinder, whereas its exterior 
counterpart will peel off. 


Whereas the second order correction to the transmitted force can have either sign (correlating with the sign of the 
cosine of the double angle), one has to proceed beyond perturbation theory to determine if the sign of the force in 
non-elliptical states changes in longer loops, as we will see in Sec. |IV B| it does. 

An interior bound figure-eight (n = 1, p = 2) would be expected to unwind into the planar Euler elastic vertical 
confined state; an exterior one, on the other hand, will constrict the cylinder. Note that this configuration possesses 
a lower energy than a free twice-covered circular loop. There is no inclination to detach from the cylinder. With 
increased length, as we see in our non-perturbative treatment of the problem, we would not expect this state to remain 
completely attached; intuitively we would nonetheless expect it to continue to constrict the cylinder. The topological 
obstruction provided by the cylinder prevents the doubly wound elastic loop from unraveling into its singly wound 
ground state. This will be confirmed when we venture beyond perturbation theory. An exterior bound loop wound 
more than once around a cylinder will always constrict the cylinder. 


2. Stability of perturbed loops 

Modulo the EL equation, the second variation of the energy is given by 

S 2 H b = j ds$£$, (56) 

where £ is a self-adjoint fourth order differential operator, defined at lowest order and about a p-fold covering of the 
equator by 


(57) 

and <f> is the normal deformation along 1, $ = <5Y ■ 1. On account of the fixed length condition, the permissible normal 
deformations ought to be orthogonal to the geodesic curvature, so $ must satisfy the global constraint 02] 

j) ds K g <f> = 0 . (58) 

The eigenmodes of C are 1, ip, cos Qip, sin Qip. s . Since the periodic eigenmodes ought to have the same number of 
windings as the original loop, the wave number is given by Q = k/p, where k is another integer. Their corresponding 


7 Its relation to the non-vanishing axial torque (mo) associated with the confined loop in this limit is somewhat mysterious. 

8 Recall that to lowest order s is given by (/?. 
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eigenvalues are 


C k = ^(tf-n*). 


(59) 


There are four zero modes with C k = 0, the constant and linear modes 1 and ip and two (cosine and sine) eigenmodes 
with the same number of periods as the original loop, k = n (Q = q). However, the eigenmode ip is not periodic and the 
eigenmode cos qip oc n g does not satisfy the global isometry condition (581, so both of these modes are not considered. 
The allowed zero modes correspond to a translation along ($t = 1) anaa rotation about (<f>^ = sing</5 oc n' g ) the axis 
of the cylinder. 


In Fig. [2] the eigenvalues C k are plotted for the first n-fold perturbations of a single covering of the equator {q = 1). 
These eigenvalues are all positive for the elliptic state with n = 1. All states with n > 2, however, exhibit negative 
eigenvalues, which signal the possibility of decay to states with lower energy, if they are accessible. Therefore, in this 
approximation only the elliptical state n = 1 is stable, whereas the higher ro-folds with n > 2 are unstable. Of course, 
to access the stability of the elliptic ground state one need to also treat perturbations lifting the loop off the surface. 
In the absence of an adhesive force one would expect this state to be unstable with respect to such perturbations. 


Ck 



FIG. 2: (Color online) Eigenvalues C k corresponding to deformations of a parallel circle. For all n > 1, the spectrum possesses 
a negative eigenvalue. 


B. Non-linear regime 

We now examine loops in their full non-linear glory. Their analysis is again facilitated by exploiting the analogy 
with a particle in a periodic potential. 


First, the constant c can be cast in terms of the turning points of U and to: let a' = 0 at a = ±a« in the quadrature 
Eq. (25 I then 


c = — - cos olm — to cos olm ■ 


(60) 

For a circular loop this reproduces Eq. ( |37[ ), relating c and to to lowest order in perturbations, c vanishes if the curve 
is vertical on the equator (% = 7 t/ 2) or when m = —1/2 cos 3 q;m (&m 7^ 7r/2). The relationship (601 permits us to 
recast the quadrature in terms of olm and to: a' 2 = V(a) 2 , where 


V{a) 2 = 2 (c — 17(a)) = (cos a — cos«m) [(cos a + cos «m) (cos 2 a + cos 2 q:m) + 2 to] . 


(61) 


Since ir > o.m > a > 0, one has cos a > cos olm > 0, so V 2 > 0 implies the lower bound on the axial torque, 
to > —2 cos 3 ctM', thus there are no closed curves with to < —2. Equation (601 now implies the bounds on c, 


c < 3/2 cos 4 OtM for 0 < ajj < n/2 (c > 3/2 cos 4 a« for 7 t/ 2 < < 7r). These bounds together imply the inequality 

c< 3/2 (|to|/2) 4 / 3 . 
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For a given p-fold covering and n-fold dihedral symmetry, Eq. (36b I establishes the relationship between m and ajf. 
Curiously, this is independent of the length of the loop. Modulo this condition, a« or m is determined as a function 
of Ar using Eq. (36a). To trace the trajectory of the loop one integrates the quadrature for a given Ar to obtain a 
as a function of l\ the relations (291 permit one to position the loop on the cylinder. 


The functional relationship between m and % is presented in Fig. [3] for each of the three sequences, L 2 p (blue curve), 
L,! (black) and L 12 (red). While this may be the first relationship to be established, the information it conveys is 
not immediately transparent. 9 It is somewhat easier to understand the behavior of m and olm as functions of A?'. In 
Figs. 0a), 0b) and gc) Q!m, m and c are plotted as functions of loop length for each of the three sequences. We 
see that m increases towards a maximum in each sequence. The corresponding loops with maximum values of m are 
illustrated in Figs. 0c), 0b) and 0d), respectively. Beyond this single maximum m decreases monotonically to zero. 
The axial torque vanishes in all long loops. In each of the L 2 p and Lip sequences, m changes sign from negative to 
positive at some finite Ar. In particular, there is a finite value of Ar at which the axial torque vanishes. This does 
not, of course, mean that the total torque vanishes. 10 These states are tractable analytically and are discussed in 
detail in Appendix [D] Surprisingly, one finds that cum does not increase monotonically with loop length. It increases 
from its value olm — 0 reaching a maximum value (> tt/2) before returning asymptotically to tt/2 in long loops. At 
some critical length on the way up, olm = tt/ 2, indicating that the loop develops vertical tangents on the equator. 
In these states the constant c vanishes, which indicates that they are also special with respect to length. Loops with 
values of olm > tt /2 develop overhangs on each side of the equator. 



a M 


FIG. 3: (Color online) Behavior of the vertical torque m as a function of am for the sequences L 2 ,i (blue curve), Lip (black 
curve) and L] .2 (red curve). The complicated relationship between m and om upon the development of vertical tangents with 
OLMi r/2 (detailed in the inset) is clarified by examining separately the functional dependencies of m and an on Ar. 


1. Lip ground state 

Equilibrium loops completing one period in one trip around the equator are represented in Fig. [5] focusing on the 
behavior as the length of the loop is increased. If Ar is small, the loop is a tilted ellipse described accurately by 
perturbation theory about the equator. See Fig. 0a). As the length increases, the eccentricity of this loop increases, 
obliging it to tilt away from the equatorial plane (Fig. 0b)); it also begins to bend out of the plane of the ellipse 


9 Note, however, that it is consistent with the perturbative behavior presented in S ect. Eva] for small aM- Recall that for a p-fold 
covering of the equator, with a = aM = 0, one has m = ( n/p ) 2 — 2. (See Eqs. ( |42| ) .) The leading correction is quadratic in aM , 
m2 = | (3 — \ (n/p) 2 ) ol 2 m . 

10 Likewise, M does not vanish in long loops. For while they are stretched along the cylinder, as we see, the asymptotic geometry is non 
trivial. 
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A r A r 

(a) (b) 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 

A r 

(c) 


FIG. 4: (Color online) (a) Behavior of (a) «m and (b) axial torque m and (c) parameter c as function of excess radius Ar for 
each of L 2 ,i (blue curve), Lyi (black curve) and Li t 2 (red curve) sequences. In all three cases a« initially increases with Ar 
from a value a m = 0, when Ar = 0. The corresponding states are loops covering the equator described by perturbation theory. 
«m reaches a maximum (> 7t/2) at a finite value of Ar. There is a critical intermediate loop length at which an = 7r/2, where 
the loop develops vertical tangents. As the length is increased olm falls asymptotically to 7r/2 corresponding to a loop stretched 
along the cylinder. The axial torque initially increases with Ar (m = ( n/p ) 2 — 2 when Ar = 0), reaching a maximum, tending 
to zero as Ar becomes large. Co = 3/2 — (n/p) 2 when Ar = 0 and c —> 0 when Ar is large. Its vanish correlates with a m = tt/2. 


(Fig- § c)) as inflections develop at the equatorial crossings (Fig. id)). This out of plane bending can be viewed 
as a consequence of the mismatch between the curvature within the two hairpins which follow the parallel circles on 
the cylinder and the increasing linearity of the sections interpolating between them. At some point the tangents on 
the equator become vertical with a« = 7r/2. With increased length, overhangs appear and the loop develops lobes 
on either side of the equator. With increased Ar, the lobes continue to grow (Fig. [f^e). 11 At some critical length 
the two lobes make self-contact at the back of the cylinder (Fig. [5jf)). Meanwhile the two segments interpolating 
between the hairpins become increasingly vertical. The loop will eventually self-intersect. Very long loops consist of 
two hairpins connected by two very long self-intersecting almost vertical sections. Both the length of these sections 
and a number of self-intersections along them is approximately proportional to the length of the loop. 


11 The angle clmi on the other hand, reaches a maximum, returning asymptotically to 7r/2. 
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FIG. 5: (Color online) Li,i loop sequence, (a) A r = 0.01 (om ~ 7r /16), (b) Ar = 0.192 («m « 7t/4 and hsm in), (c) A r = 0.78 
(cum ~ 27r/5 and m m ax), (d) Ar = 1.8 (q m = 7r/2 and /ismax), (e) Ar = 5.544 (cum max ~ 0.548 n) and (f) Ar = 27.819 
(om « 0.522 7t, first self-contact occurs). In (f) only the upper half of the loop is displayed. The normalized magnitude of the 
normal force A is color coded along the loop in these figures. 


2. Li ,2 excited states 

Significant features of loops undergoing two periods while they wrap the cylinder once are represented in Fig. [6] 
If the loop is short, it consists of a small oscillation about the equator, (Fig. [6](a)), the amplitude of the oscillation 
increasing with length (Fig. [6](b) and (c)). When Ar = 1.313 the loop becomes vertical at the equator, (Fig. [6](d)). 
Lobes develop in longer loops, with olm reaching its maximum when Ar = 5.48 (Fig. |6je)). The first self-contact 
occurs when Ar = 8.59 (Fig. if))- Self-intersections will occur in longer loops. Asymptotically, states consist of four 
hairpin bends interpolated by increasingly vertical mutually intersecting sections. One can also think of this state 
as a folded n = 1 state, with the folds themselves forming hairpins. In general, states with n > p will behave in a 
manner similar to those in this sequence. 


3. L2 ,i states 

Non-trivial loops which wind about the cylinder more than once (p > 1) are also possible whenever r > p. The 
L- 2 a loop sequence is represented in Fig. [7J If the excess length is small the loop approximates a double covering of 
the equator, adopting as the length is increased the shape of a folded figure of eight wrapping the cylinder, with its 
self-intersection lying on the equator. See Fig. [7^a) . (If p = 3, the two intersections migrate off the equator.) The 
crossing angle increases as the figure is stretched. See Figs. [7](b)-(d). This angle becomes vertical when A?’ = 11.06, 
see Fig. [7je). olm reaches a maximum value 0.51 n, at Ar = 27.548. See Fig. [7]T). In longer loops, additional 
self-intersections appears above and below the equator. In contrast to the asymptotic behavior of the \ sequence, 
the two hairpins lie on the same meridian. 

The appearance of the self-intersections in these three sequences of loops can be understood more easily by unfolding 
them as described in detail in Appendix |E) 
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FIG. 6: (Color online) Li^ loops sequence, (a) A r = 0.01 («m ~ 7r/16), (b) Ar = 0.13 («m ~ 0.216?r and m ma x), (c) 
A r = 0.365 (om = 7t/ 3), (d) Ar = 1.313 (a« = tt/ 2 and hs max), (e) Ar = 5.48 (aitmax ~ 0.595 7r) and (f) Ar = 8.59 ( first 
self-contact occurs, ^ 0.588 n). The normalized magnitude of the normal force A is color coded in these figures. 


4- Trajectories in the c-m parameter space 

For completeness, we also represent the three sequences as trajectories in the c — m parameter space. These 
trajectories are bounded by two curves, 

1. The line m = —1/2 — c, denoted by II, representing loops covering the equator p times, is indicated by the gray 
dashed line in Fig. [8j 

2. The curve c = 3/2(|m|/2) 4 / 3 , denoted T representing loops saturating the inequality between c and m, is 
indicated by the gray dotted curve in Fig. [8] 

The point (3/2,—2) at the intersection of II and T, is a limit point, corresponding to loops winding the equator an 
ever increasing number of times (p —¥ oo), with a finite number of periods, so that cq —¥ 3/2 and mo —> — 2. 

All three trajectories originate on II and terminate at the origin. Infinitely long loops have vanishing m and c. In 
this limit the fixed length constraint gets relaxed as c becomes vanishingly small. 
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FIG. 7: (Color online) Loops of sequence L 2 ,i. (a) A r = 0.02 (oim ~ 7r/16), (b) Ar = 0.25 (om ~ 0.1947r), (c) Ar = 1.149 
(um = 7t/3), (d) Ar = 6.916 (oim ~ 0.487T and m max ), (e) Ar = 11.06 (om = 7t/2 and /ismax) and (f) Ar = 27.548 
(«m max ~ 0.517r). In (f) only the upper half is displayed. The normalized magnitude of the normal force A is color coded in 
these figures. 


The trajectory for Lip (I^i) indicated by the black (blue) curve in Fig. [8]begins at the point (1/2, —1) ((5/4, —7/4)) 
spiraling towards the origin. Both trajectories cross the line c = 0 three times; twice with a finite value of m, one 
with an acute angle and another with vertical tangents at the equator, and again at the origin. 

Likewise, the trajectory for Lip begins on II at (—5/2, 2) but it forms an arc, crossing the line c = 0 twice, once where 
the loop has vertical tangents and again at the origin. 

In general, any sequence with n > p begins on II to the left of Li p, whereas sequences with n < p begin on the right, 
all sequences tend asymptotically to the origin. Sequences with c > 0 (c < 0) or n/p < y/3/2 ( n/p > 3/2) will cross 

the line c = 0 three times (twice) and thus will have spiral-like (arc-like) trajectories as in sequence L2 ,i (Lip). 


5. Total energy 

Loops in the L Pjra sequence originate in p-fold coverings of the equator. Thus their initial energy is hs = 7r p, 
independent of n. See Fig. [9] If the excess radius Ar is small, loop states in the Lip sequence are approximately 
elliptic; their energy is thus close to that of their planar vertical counterparts presented in Appendix [C] as revealed by 
the initial coincidence of the solid black and gray curves, respectively, in Fig. [9] However, in this sequence the energy 
reaches a global minimum when Ar = 0.192, corresponding to the state illustrated in Fig. (5Tb) , increasing towards a 
global maximum value hs = 1.042 tt when Ar = 1.80 (the state with o?m = 7t/ 2 and c = 0 J, after which it decreases 
asymptotically towards Hb = 7r. Long loops and short ones have the same energy! By contrast, the energy of planar 
vertical loops decreases monotonically to a minimum value of 0.914-7T when Ar = 0.393 which is when the curvature 
in the free loop vanishes where it makes contact with the wall. Thereafter it remains constant because any additional 
length gets directed into the straight line segments. See the inset in Fig. [9j We see that vertical loops always have 
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C 

FIG. 8: (Color online) Trajectories in the parameter space c—m of sequences L/ 2,1 (blue curve), Lip (black curve) and Li^ 
(red curve). Dashed line II represents the bound m = —1/2 — c, corresponding to a single or multiple covering of the equator, 
whereas dotted line T represents the bound c = 3/2(|m|/2) 4 ' 3 and correspond to loops covering the equator a very large number 
of times. The intersection of II and T occurs at the point (3/2, —2). The inset shows the trajectories in the neighborhood of 
the origin, corresponding to loops with very large excess radius. 


lower energy than any bound counterpart: the former thus provides the ground state of interior bound loops allowed 
to unbind. 

The energy of excited states in the Li i2 sequence is indicated by the red line in Fig. [9j While it is initially degenerate 
with the Liq sequence, with energy hs = 7r, unlike the latter it increases with the excess radius from the beginning, 
reaching a maximum at Ar = 1.313 (also when olm = 7r/2), thereafter decreasing asymptotically to the value hs = 2 n. 
The energy of the double winded L 2j i sequence is indicated by the blue curve in Fig. [9] It starts at hs = 2 7r; unlike 
the energy of the L^i sequence which only decreases initially, the energy decreases monotonically to hs = 7r, the 
same asymptotic energy as a loop in the Lip sequence. 

In general, for L p>n sequences with n > p, the maximum of Hb occurs in states with vertical tangents on the equator 
and finite vertical torque, i.e., ajj = 7r/2 and m ^ 0. The corresponding states in the Lip and Li j2 sequences are 
displayed in Figs. id) and [5 e). For n < p the initial states possess the maximum energy within the sequence. Since 
loops with very large excess radius are essentially straight almost everywhere, the principal contributions to the total 
energy will be from the 2 n curved hairpins, each with approximately the same energy as a semicircle, Hb = 7r/2. This 
explains the coincidental simplicity of the asymptotic energy in the three sequences: for all L PiTl , the energy hs —>• 7rn 
(from above) as Ar —► 00 . It is independent of the topology of the loop. 


6. Transmitted normal force 


The magnitude of the normal force A, scaled to lie between —1 and 1, is represented in color along the loops in 
Figs. [5j[6]and[7] indicating clearly the regions where they are pushing or pulling the cylinder: red (purple) represents 
regions where the push (pull) is maximum, 12 whereas green represent regions in which the transmitted force is low. In 
general, it is observed that the maximum push occurs at the equatorial crossings, whereas the minimum (maximum 
pull if negative) occurs at the hairpins. The maximum (solid curves) and minimum (dashed curves) values of A are 
plotted as functions of the excess radius Ar in Fig. 0a) for each of the three L^i (black curve), L 1?2 (red curve) 


12 


Recall that this convention is adapted to the viewpoint of an interior loop; if the loops is outside, push and pull are reversed. 
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FIG. 9: (Color online) Total energy of loops belonging to the three sequences L, 2 ,i (blue), Lip (black) and Li i2 (red), as well 
as the energy of the vertical loop touching the cylinder along its sides (gray). The initial value of the total energy of the L P:Tl 
sequence is hs = 7r p, corresponding to a p-fold covering of the equator, whereas the energy of loops with very large excess 
radius tend asymptotically to the value hs = 7rn, indicated with dashed lines. 


and L 2j i (blue curve) sequences. 


The perturbative result, Eq. (551, indicates that the initial state of a sequence L p rl exerts a normal force of magnitude 

A 0 = ( n/p) 2 - 1. (62) 

Regardless of n and p, the maximum and mini mu m values of A tend to A max —> 1 approximately and A m i n —> — 2 
exactly as Ar —> oo. In fact, in this limit Eq. (341 reads A « 2 cos 4 a (5 — 6cos 2 a), with minimum at a = 0 and a 
maximum when cos 2 a = 5/9, so that X ma x = 250/243, which is 1 for all practical purposes. 


For completeness, we present the total transmitted force A = £ ds A, which indicates, whether, on average the loop 
is pushing or pulling on the cylinder. It is plotted in Fig. [l0|(b) in each of the three sequences. Since the p-fold 
coverings of the equator exert a constant normal force, the initial total force is simply Ao = 27rpAo- It is found that 
A tends asymptotically to Aoo — > n n from below as A?’ — > oo. Whereas enclosed bound loops both push or pull on 
the cylinder, depending on the sign of Ao, ultimately all such loops end up pushing the cylinder. 


In the Lip sequence, the single covering of the equator does not initially exert any force, Ao = 0. As Ar increases, 
however, second order corrections kick in with A alternating from positive to negative along the loop. This behavior 
persists with increasing Ar; see Fig. [5] and the black curves in Fig. [lOja). The negative pulling force on an inner 
bound loop is increasingly localized in the neighborhood of the tips. We saw earlier that the energy of this state is 
always higher than that of the planar vertical loop of equal length spanning the interior. Importantly, however, the 
averaged force is always positive. While the energy would suggest that the Lip ground state is always unstable, the 
sign of the contact force suggests otherwise. What appears to occur is that the loop unbinds from the cylinder in 
the neighborhoods of the tips. If it is short it will rotate into the vertical. As the length increases, however, the 
lobes pushing on the cylinder snag the loop providing a potential barrier obstructing the passage towards the planar 
vertical loop. A detailed perturbative study, beyond the scope of this paper, is required to settle the issue. 


Along loops in the Li „ sequences, initially A is uniformly positive (see Fig. 6[a)), indicating that an Euler type 
instability is associated with their formation as interior bound states. However this pus h on the cylinder persists only 

7.65xl0 -2 for n = 1 (see the red dashed line in Fig. 10 a)), the minimum of A turns 


for short loops; once Ar > 7.65x10 2 for n = 1 (see the red dashed line in Fig. 
negative within the regions near the tips, indicated in Fig. §b). In the case of an interior bound state, one would 
expect the loop to unbind in the neighborhood of its tips so as avoids these regions. 


An exterior bound loop with p = 1, free to unbind will completely unbind, reforming as a circular loop. As we have 
already glimpsed, loops in the L 2j i sequence behave very differently: now A is initially uniformly negative. While an 
inner bound loop in this state is manifestly unstable initially, its exterior bound cou nterpart is stable. However when 
Ar > 0.317, the maximum of A becomes positive (c.f. the blue solid line in Fig. 10 a)) in the regions surrounding the 
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FIG. 10: (Color online) (a) Maximum (solid curves) and minimum (dashed curves) of the magnitude of the normal force A for 
sequences L 2.1 (blue curve), Li^i (black curve) and Li _2 (red curve), (b) Total normal force of the same sequences. 


center of the hairpins; see Fig. [7] Beyond some higher value of A?’ the average force becomes positive. This suggests 
that a sufficiently long interior loop may be stable modulo shortcuts. 

As discussed in the perturbative context, doubly-wound exterior loops are more interesting because such loops will 
always constrict the cylinder. If the loop is short, it will constrict the cylinder everywhere, albeit not in an axially 
symmetric way. If longer, the contact force will favor the unbinding of the exterior bound loop in the neighborhood of 
the crossings. However, such loops will remain attached to the cylinder at its hairpins applying a force on the cylinder 
at the points of contact. In any case the cylinder provides a topological obstruction preventing the doubly wound 
elastic loop from unraveling into a singly-wound circular loop. This happy conspiracy permitted by the topology has 
no analogue on a sphere. 


V. DISCUSSION AND CONCLUSIONS 

We have examined the equilibrium shapes adopted by an elastic loop either confined by a cylinder or constricting 
it. Our initial motivation was to contrast the behavior of loops confined within a cylinder with that of their spherical 
counterparts, and specifically to explore how these loops respond to the changed topology. Loops may wind a cylinder 
any integer number of times. They may oscillate one or more times as they do this. Whereas the winding has no 
topological significance for a loop that is free to unbind within the cylinder, it very much does if it is bound to the 
outer surface. 

On a sphere, the ground state of a confined loop with a radius exceeding the spherical radius is attached to the sphere. 
No additional adhesive forces are necessary. The loop pushes against the sphere and it will remain bound everywhere 
even if free to detach. Instabilities of excited states occur through rearrangements on the sphere itself. When there 
are topological obstructions, they may be bypassed by stepping into the interior. Elastic loops, however, simply do 
not bind to the exterior surface of a sphere. What occurs on a cylinder is a lot more interesting. 

On a cylinder, one needs to distinguish between stability on the cylinder and off it. We have considered, in detail, the 
behavior of states bound to the cylinder. If the loop is short there is a single bound ground state for each winding 
p, the state with n = 1. All states with n > 2 will be unstable if the loop is short. However, as loops grow in length 
the sections connecting the troughs and crests straighten out along the axis; if the loop oscillates more than once, the 
oscillations morph into folded hairpins. All loops grow by extending indefinitely in the axial direction, which permits 
them to relax their axial torque. One can show that for each n > 2 there will be a critical length above which the 
bound folded loop is stabilized. There will be a crossover from instability to stability associated with the exploration 
of the long (or non-compact) direction. This behavior has no analog on spheres. 

If states are free to unbind, the behavior inside and out are very different; both are interesting. Unlike loops binding 
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to spheres, the contact forces are not positive everywhere. Indeed, they are always negative in the neighborhood of 
the crests and troughs of the oscillations. Interior bound states would be expected to unbind within these regions, 
lowering the loop curvature and, as a consequence, lowering the loop energy. The ground state itself with p = 1, 
n = 1 will also unbind, and if sufficiently short, will be unstable with respect to deformations rotating it off the 
surface into the confined ground state: a planar loop compressed along its sides. Experimenting with a loop in a 
cylindrical wastepaper basket or pipe confirms this prediction. Longer partially bound states may become stable 
beyond some critical length. Folded hairpins will certainly be stable modulo minor adjustments. A stability analysis 
accommodating normal perturbations is necessary to settle these issues, but it is not going to be simple. 

We saw that states with p > 1 and n = 1 would tend to unbind if below some critical length. This is because the 
contact forces are inward everywhere. As these states get longer the force turns to push except at the tips so that 
the state would be expected to stabilize modulo the same sort of corner-cutting we have described. These states are, 
however, far more interesting placed outside the cylinder. 

An exterior loop wrapping a cylinder once will always unbind in the absence of an adhesive force holding it there, 
reconfiguring itself as a simple circular loop. However, one that wraps the cylinder more than once is constrained 
by the topology to remain in contact with the cylinder, completely initially and partially-in the neighborhood of 
its tips-beyond some critical length. In either case it will always apply a compressive force on the cylinder in these 
neighborhoods. This force generally is not axially symmetric. This behavior can also be confirmed in a broom 
cupboard experiment by coiling an elastic loop into a figure-eight and wrapping a broom handle: whereas a short 
loop will hug the handle completely, a long one will detach everywhere except at four points. This experiment also 
reveals that the detached state does not possess the symmetry of its bound counterpart. Indeed, if the loop is coiled 
three or more times, it is clear that the sequence of crossing (over or below) is not unique and the detached state 
will depend on the particular sequence. This mechanism may be relevant to an understanding of the constriction of 
membrane necks in dynamin mediated membrane vesiculation. It should be remarked that we never anticipated the 
possibility that loops wrapping a cylinder could bind without any additional agency, never mind, finding that they 
provide a possible mechanism for neck constriction. 

A membrane neck is, of course, better modeled as a catenoid. Modeling the constriction of a catenoid is technically 
more challenging because, along with the surrender of the translational symmetry along the cylinder along its axis, 
one also cedes the quadrature. 

An obvious direction for future work is to study open elastic spirals on cylinders, catenoids and other axially symmetric 
geometries with free boundaries. Spiral geometries represent the behavior of chains of assembled proteins such as 
dynamin. However, unless physically unrealistic BCs are introduced, the simple elastic energy considered here will 
not constrict the geometry. A constricting spiral will require an environmental bias, such as a spontaneous normal 
curvature, if it is to apply a force. If it is to spiral, it will require a constraint or bias on the twist which translates 
into a constraint or bias on the geodesic torsion. This still leaves us with a long way to go before we are in a position 
to treat a spiral binding to a fluid membrane taking their mutual interaction into account. The consolation is that 
there is interesting physics to be picked up along the way. 
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Appendix A: Axisymmetric Surfaces 

To describe an axisymmetric surface, it is convenient to introduce cylindrical coordinates ( R , <p, Z), adapted to the 
symmetry with tangent vectors 


p = (cosip, sinip, 0) 


(p = (— sin ip, cos p, 0) , 


z = (0,0,1). 


(Al) 
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The surface is now parametrized by arc-length l along a meridian and the azimuthal angle ip along the parallel, 

x(e,<p) = R{e)p + z(e)z, (A2) 

where R(£) 2 + Z(£) 2 = 1 and the dot represents derivation with respect to l. The two tangent vectors adapted to this 
parametrization are 


e e = Rp + Zz, e v = Rip . (A3) 

Thus, the line element on the surface is given by ds 2 = <M 2 + R{t) 2 dip 2 , and the metric tensor assumes the form 

9ab = ( 0 R 2 ) ' 

The outward unit normal vector, n = <? -1 / 2 x e^, is 

n = Zp — Rz ; (A5) 

the extrinsic curvature tensor is diagonal, 

A '-‘ = ( ~T RZ ) ' < A6 > 

The eigenvectors of the shape operator K a b = g ac K c t (its principal directions) lie along meridians and parallels, with 
corresponding eigenvalues given by 

kj_ = —R/Z, k\\ = Z/R. (A7) 

These are the curvatures of these curves. The two symmetric invariants of the shape operator are the mean and 
Gaussian curvatures K = + k || and Kq = k± kh. 


Appendix B: Hamiltonian approach to cylindrical confinement 


Here we provide a direct derivation of the “second” integral of the EL equation describing the confinement of a 
semi-flexible polymer within a cylinder. The “Lagrangian” density 2 z? is constructed in terms of the three coordinates 
[ip, z, a ) 13 characterizing the curve T. These variables are not independent; they are related by Eqs. (12 1 . It will be 
necessary to implement these constraints using two local Lagrange multipliers X v and A z . One thus constructs the 
following effective energy 


Jf[a,ip,z,a',ip',z',X ( p,X z \ = j ds + — 7 j— + ~ cos a) + X z (z' - sina)^ . 


The momentum densities conjugate to the coordinates (ip,z,a), are given by 

Pz = A z , P a = OL . 


Pip Xip, 


(Bl) 


(B2) 


Since <p and 3 are cyclic coordinates (because of the rotational and translational symmetries) P v and P z are conserved. 
Also, L does not depend explicitly on arc-length. Thus, its Legendre transformation R = ip 1 P^ + z'P z + a'P a — C is 
constant. Evaluating it, one obtains 


1 2 cos 4 a 


+ P z sin a + Pu, cos a . 


(B3) 


This reproduces the “second” integral Eq. (24) with the identifications, H = c, / = —P z and m = — P^. 


13 As before, here we scale lengths with the cylinder radius. 








Hamilton’s equations for a are 
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dP 

dP n 


= P n 


P' = 


dT~L 

— —— = —2 cos 3 a sin a + P v sin a — P z cos a . 
da 


(B4) 


To reconstruct the curve one solves these equations with initial conditions a(0) = ao, -P Q (0) = P a o and fixed P v and 
P z . The azimuthal and height functions are obtained from integration of the two Hamilton equations 


7 

ifi = 


dP 

dPv 


cos a , z' 


dP 

dP z 


sin a. 


(B5) 


Without loss of generality one can use the initial conditions </?(0) = 0 and z(0) = 0. The three constants P z , P v and 
P are determined from boundary or periodicity conditions. 


Appendix C: Confined planar curve 

Let us consider the confinement of a closed curve of radius R, lying vertically inside a cylinder of radius i?o, touching 
the cylinder tangentially on its two sides. We have claimed that this state has lower energy than its bound counterpart. 
In the absence of additional adhesive forces, the bound interior ground state-if sufficiently short-would be expected 
to be unstable with respect to deformation into the cylinder rotating it into this state. Since the curve is planar, it is 
described by the EL equation mm 


£n = k" + k — cj = 0. 

The translational invariance of the bending energy permits one to identify a quadrature 



(Cl) 


(C2) 


where F is the magnitude of the force vector F along the planar curve, given by 

F = ( - c ) T + k'N . 


(C3) 


F is conserved along the curve, i.e., F' = 0. The quadrature (C2) is integrable in terms of elliptic functions maun. 
The specific function will depend on the relative values of F and c |TB] [BO] ■ 


If F 2 < c 2 the FS curvature is given by 

k(s) = 2qdn[qs,u] , F = uq 2 , c = q 2 (2 — u), (C4) 

where the Jacobi elliptic function with argument cf> and parameter u is defined by dn 2 [(f>,u} = 1 — wsn 2 [<j>,u] and 
sn[<j>, u] is the Jacobi elliptic sine |45j . In one period the curvature will oscillate asymmetrically about the value \/2c 
between two positive (or negative) values, so one has 0 < n m \ n < k < K max , where K m in = 2 q \J 1 — u and K max = 2 q 
(0<u<l, so F > 0 and c > 0). Thus the curve has an orbitlike shape. 


If F 2 = c 2 the curvature is given by 


k(s) = 2gsech qs , F = c = q 2 . (C5) 

The curve is not periodic; the curvature is positive everywhere but it takes an infinite length for the curve to become 
planar. That is k —> 0 as s —> ±oo. This curve is termed borderline, corresponding as it does to the separatrix 
separating orbital and wavelike behavior. 


If F 2 > c 2 the curvature is given by 

k(s) = 2 y/u q cn[< 7 S, u] , F = q 2 , 


c = q 2 (2u- 1), 


(C6) 
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where ci\{(b, u\ is the Jacobi elliptic cosine [33]. The curvature oscillates symmetrically about 0 between ±Km, he., 
—km < k < Km, where km = 2 y/uq (0 < u < 1), so the curve adopts a wavelike shape. 

The parameters q and u are determined by specifying BCs, which requires knowledge of the embedding functions of 
the curve. 

The Cartesian coordinates of the curve Y = (x(s),y(s)) and its tangent vector T = (x'(s), y’(s)) can be expressed in 
terms of the curvature. Since F is a constant vector, it is more convenient to work in Cartesian coordinates and align 
it along direction x, soF=Fx, with F constant. Projecting F onto the Cartesian basis (x, y), yields 

x' - cj - Z/V = F, y' + x 'k' = 0. (C7) 

Taking linear combinations of these equations, we get the components of the tangent vector in terms of the curvature 
and its derivative, 


. ru 

Fx = — c, 


F y' = —k! . 


(C8) 


This can be regarded as a factorization of the quadrature (C21. Introducing the angle ip that the curve makes with 
the direction of x, the components of the tangent vector can be written as x' = cosip and y' = simp. Also, the 
curvature is given by the derivative of this angle, k = ip'. Equations (C8l can now be mapped into the statement 


of conservation of energy E and the equation of motion of a pendulum: arclength is identified with time and ip is 
identified with the angle that the pendulum makes with the vertical (see, for example, USD- 

Let arclength be measured from the midpoint of the curve, so x(0) = 0. Let the contact points occur at an arclength 
s = is},. The normalized loop length is l = 4 S/, = 2irr, where r = R/Rq. Confinement places two constraints on the 
curve. First, since the curve must align with the container at points of contact, the tangent vector must be vertical 
at liftoff, therefore x'(±Sb) = 0 and y'(±Sb) = Tb Thus, from Eq. (C8) one has that the constants c and F are given 
by the curvature and its derivative at the contact points, 


c(±Sb) = y/2 c (c > 0), k' (±Sh) = zk F . 


(C9) 


The first condition permits one to obtain the angular wavenumber q in terms of the modulus to. Using expressions 
(C41, (C5) and (C6l of k and c in terms of q and u, this condition read for each case 


L4’ U 


qsb = F 
qsb = arcsinh 1, 
qs b = F 


arcsm 


yplu 


F 2 < c 2 , 
F 2 = c 2 , 
F 2 > c 2 . 


(ClOa) 

(ClOb) 

(ClOc) 


where F[</>, u] is the incomplete elliptic integral of the first kind, |43j . This condition is of an intrinsic character, for it 
determines the distance of the contact points along the curve. The second condition due to the confinement on the 
curve is of extrinsic character, for it fixes the spatial distance between the contact points: the coordinates x of the 
contact points are constrained to be equal to the cylinder radius, 


:-(±Sb) = ±Rq 


(Cll) 


To implement this constraint, one needs to integrate Eqs. (C8l. The second equation in (C8l can be integrated 


readily, which is a consequence of the rotational invariance about the direction orthogonal to the plane of the curve. 
Indeed, the integration constant is identified as the conserved torque along direction z; that is, Fy = —k — M, where 
M = M • z. By setting the contact points at y = 0, it follows that the magnitude of the torque is given by the negative 
of the curvature at the contact points, M = —Kb■ Therefore, the coordinates of the curve are given by 


S 

x(s) = J dt(^ 


K(t) 2 


— C 


2 / 0 ) = ^ Ob - k{s)) . 


(C12) 
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In full, the Cartesian coordinates for each case read 


x = — Efamfgs, u], it] + ( 1- 

uq ' 

2 

x = - tanh qs — s , 

q 


x = - E [am[qs, zz], it] — s , 


2 

u q 
2 


y= — ( A/ 1 — dn[gs,u] 


V = 


y = 


1 

q V \/2 
2 


— sech qs I , 


\J U ~ \ ~ cn[qs, zz] j 


F 2 < c 2 ; 

(C13a) 

F 2 = c 2 ; 

(C13b) 

F 2 > c 2 . 

(C13c) 


where E[</, zz] is the incomplete elliptic integral of the second kind and amj^a] is the Jacobi amplitude [53]. This 
segment forms the bottom half of the free loop; the upper half is obtained by an appropriate up-down reflection. 


Imposing condition (Cll I in Eqs. (C13a) - (C13cl provides another expressions of the wavenumber q in terms of the 
parameter m, obtaining for each case 


Rq u q = 2 E 


L4 ,U 


- (2 — zz) F 


L4’ U 


i?o q = v^2 — arcsinh 1, 


Rq q = 2 E 


arcsin 


\/2 zz 


,u 


-F 


arcsin 




,u 


F 2 < c 2 , 
F 2 = c 2 , 
F 2 > c 2 . 


(C14a) 

(C14b) 

(C14c) 


By combining Eqs. (CIOI and (C14l and using the relation s& = zrr/2 one determines the excess radius A?’ = r — 1 
as a function of u 


_ uF[j ,u] _ 

^ (E[f,«] - (1-f) F[f,«]) 

2 arcsinh 1 


7r (-s/2 — arcsinh l) 


F 

arcsin /—, u 

\J2u 


7T 

( E 

arcsin — h=, u 

V 2ii 

- -F 
2 l 

arcsin — h= , u 

\/2 u 

) 


F 2 < c 2 , 

(Cl 5a) 

F 2 =c 2 , 

(C15b) 

F 2 > c 2 . 

(C15c) 


This makes it possible to parametrize the loop by u. Alternatively one can solve Eqs. (C15) numerically in each case 
to determine u for a given excess radius Ar. 


The total energy of the loop is easily obtained by noticing from Eq. 
x' . Thus, the energy of one quarter of the loop is given by 


(C81 that the energy density is proportional to 


Hb 

4 



Sb 

j ds (F x' + c) = F Rq + c Sb , 
0 


(C16) 


which establishes the role of the constants F and c as the conjugate variables to the extrinsic and intrinsic lengths 
of the curve respectively (cf. Eq. (32) for the confined case). Using Eqs. (ClOl and (C14), one gets that the total 
energy of the loop in each of the three cases by 


H b = 
H b = 


8?E L4’ 
4^2 q, 


H B =8q 



1 

arcsm ,_ , u 

V2u 


(1 — u) F 


arcsin 


1 


x/2 u 



F 2 < c 2 , 

(C17a) 

F 2 =c 2 , 

(Cl 7b) 

F 2 > c 2 . 

(C17c) 


One quarter of the confined loops are plotted for increasing radius in Fig. mt a). If u = 0 one has a circle. In the range 
0 < A?’ < 0.053 (0 < zt < 1) the curve is consists of orbitlike elastica segments, whereas for 7 > 0.053 (1/2 < u < 1, u 
decreases as the length increases) it consists of wavelike elastica segments (at Ar = 0.053 (u = 1) consists of borderline 
elastica segments). In this wavelike regime, when Ar = 0.393 (u = 1/2) the curve adopts the limit shape shown in 
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Fig. |ll[ a) and by further increasing the length, the middle region of the loop-where the curve makes contact with 
the cylinder-will elongate, so the two limit arcs will be separated by straight lines. The energy of the loops is plotted 
in Fig. E3b) as a function of their radius. The bending energy Hg of the confined loops begins at Hg/ir = 1 /Rq 
(circle with u = 0) and decreases monotonically as the radius of the loop increases, saturating when the limit shape 
with straight segments is reached for such lines do not contribute to the bending energy. 14 The normal force acting 





(a) 


(b) 


(c) 


FIG. 11: (a) One quarter of vertical elastic loops of different radii R confined inside a cylinder of Ro (Ar = R/Ro — 1). The 

loops have vertical tangent at two points where they touch the cylinder. The shortest curve is the circle with radius equal 
to that of the cylinder (Ar = 0). In the interval 0 < Ar < 0.053 the loop is composed of arcs of orbitlike elasticas (light 
gray curve), when Ar = 0.053 of a segment of the borderline elastica (dashed curve) and after which is formed by wavelike 
elastica segments (dark gray curves), reaching a limit shape at Ar = 0.393 (black curve), such that beyond this point the 
excess length will be located at the contact regions, so the vertical intermediate region will increase, (b) Bending energy of 
the confined loops as a function of their reduced radius. The black dot represent the energy of the loop formed by borderline 
elastica segments; curves with smaller (greater) radius consist of orbitlike (wavelike) elasticas. Hb starts at a value of tt/Rq 
and decreases monotonically (black curve)-but much slower than the bending energy of a free loop (gray line)- as the loop 
radius increases, reaching a limit value of 0.914 n for the curve with Ar = 0.393, after which it remains constant, regardless of 
the loop length, (c) Magnitude of the normal force transmitted to the cylinder at the points where contact is made. 


on the point where contact with the cylinder occurs is given by A = 2 F (each arc contributes with an amount F). 
This can be seen by considering a stress source of magnitude A acting on the point of contact at s = sj,, in the EL 
equation (Cl I, 


£N = k" + ft ( —-c ) = A S(s — Sb) ■ 


(C18) 


Integrating this equation in the neighborhood of the contact point and taking into account that the curvature is 
continuous across the contact region, 15 one has 


/ 


s b +Ss 


ds £n = k'i — k'_ = A, 


Sb~5s 


(C19) 


where k', and k'_ are the limits of the arclength derivative of the curvature, on the top and bottom sides of the 
contact point. Thus, the force exerted on the loop is reflected by the discontinuity in the derivative of its curvature. 
This is analogous to the force exerted on a fluid vesicle by a ring, presented in Ref. | 


From the BCs. (C91, one has 


14 This decay is slower than that of the energy in a free loop Hgo = rr/ft, and achieves a limit value rather than tending to zero for very 
long loops. 

15 If an adhesion energy is considered, the curvature would have a discontinuity 031H3- 
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k' ± = ±F, so A = 2 F. Moreover, from Eqs. (C4 (-(C6) follows that A is proportional to the squared wavenumber, 
which is plotted in Fig. |TT|c). As the loop radius increases, A increases from 0 (circle) to a constant value 1.435 (limit 
arc). When the straight lines appear, they do not exert a force on the cylinder and the corresponding normal force 
for each arc A = F is non vanishing only at the points of detachment. 


Appendix D: Closed curves with zero axial force and torque 

Loops of special interest are those in which not only the axial force vanishes but, in addition, the axial torque 
vanishes, i.e., m = 0. 1 ' While this condition is met in infinitely long loops, rather surprisingly it is also met non- 
trivially for loops of specific length and precise values of n and p. Now the constant c given in Eq. (601 simplifies to 


the scaled normal curvature as 


c = — cosr M /2 = KT n M /2, so the quadrature (241 reduces to a' 2 = cos 4 a — cos 4 ajf, which can be recast in terms of 


R' =2 


yjn n (1 - K n ) (n 2 n - K 2 nM ) . 


(Dl) 


Solutions exist in the range 0 < R u m < R n < 1, so that the normal curvature will oscillate between a minimum 
RnM € (0,1) at points where the loop crosses the equator (a = olm) and 1 at points with minimum and maximum 
height (a = 0). One finds an exact analytic solution upon integrating Eq. (Dl I, 


1 w 

arcsin ^ 

I 1 K n ^K n M 

1 KnM 

V 2 K nM 

/ 1 K n M K n + /^nM 

2 


(D2) 

Here arc-length s is measured from the tips where n n = 1. The arc-length measured from tip to crossing point is 
given by sm = K[(l — R u m)/ 2] / x/2 R n M, where K[it] is the complete elliptic integral of the first kind [53], so that the 
total length of the loop is l = 4ns m and the excess radius is Ar = 2ns m/k — P- 


The normal curvature as a function of arc-length is identified by inverting 

_ , , RnM dn 2 U/2 R nM S, (1 - R nM )/ 2] 

K n (SJ = --- 2 f _ - -----y . 

1 + K nM - dn“ v' 2 KnM S, (1 - K nM )/2\ 


(D3) 


Both the height and azimuthal functions can be determined by combining the quadrature (Dl) with the relations 
p' = cos a = and z ' = sin a = \/l — R n and integrating: 



arcsm 


-F 



The maximum height z max occurs at the tips where K n = 1. The total azimuthal extension of the loops is 


A p = 


4n 




:K 


K"nM 


1 — Kr 


M 


1 


' ^nM 


Closure A ip = 2np implies that the corresponding wave number is given by 

H /. :——- j-r r 1 n n M 

q = - = 7T vi + Kwm/k ——— 

P 2 L 1 + K nM 


(D4a) 

(D4b) 

(D5) 

(D6) 


In Fig. 12 q is plotted as a function of KnM• We see that the wavenumber is bounded in the interval 0 < q < y/2. 


Excited states with vanishing torque do not exist in any sequence with p = 1; the only such state is an elliptic ground 


16 The change of sign between the up and down regions stems from the fact that it is obtained by a reflection of the top arc, rather than 
a periodic continuation of the corresponding elastica. 

17 Note, of course, that the total non-conserved, torque does not vanish in these states. 
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FIG. 12: Allowed wavenumbers for loops with zero axial torque and force. The loops are defined for wavenumbers given by 
rational numbers in the range 0 < n/p < \/2. 



(a) (b) (c) 


FIG. 13: (Color online) Curves with zero axial force and torque with (a) p = 3, n = 4 and Ar = 0.122 (cum = 0.1257r); (b) 
p = 1, n = 1 and A r = 0.297 («m = 0.2967t); (c) p = 2, n = 1 and Ar = 4.8 (om = 0.46l7r). The normalized magnitude of the 
normal force A is color-coded in these figures. 


state with n = 1 and m = 0 (cf Fig. [3|. For a given q one has to solve numerically Eq. (D 61 to determine R u m- 
Loops with q = 4/3,1 and 1/2 are illustrated in Fig. 13 Their corresponding scaled normal curvatures at the crossing 
points are R u m = 0.853, 0.359 and 0.015 respectively. 


Loops with R n M ~ 1 will be close to p-coverings of the equator. As R n M decreases towards zero the loops elevate, 
getting more and more vertical as R u m tends to zero. In the limit case R u m = 0, the reduced normal curvature and 
the coordinate functions assume a simple form. For c = 0, the quadrature reduces to a' = cos 2 a, so tana = s and 
Rn = l/(s 2 + 1), whereas the azimuthal and height coordinates for this case can be represented by ip = arcsinhs 
and z = \Js 2 m + 1 — Vs 2 + 1, which together permits one to express the normal curvature and the height in terms of 
azimuthal angle as R n = sech 2 p and z — cosh ipM — coshp. Thus it takes an infinite number of windings and infinite 
height to go from the tips with R n = 1 (<p = 0) to the points crossing the equator with R n = 0 (<p —► oo), where 
the loops have vertical tangents. Therefore the loop has infinite excess radius and winds infinitely many times while 
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completing one period, so n = 1 and p —» oo, and consequently q = 0. 

The energy of a curve with length l = 2 Sb is given by 

h B = ^n J dsR 2 n -R 2 nM l/2. (D7) 

o 

There is no analytic expression in terms of elliptic functions, so it has to be integrated numerically. The energy of the 
limit loop with k u m = 0 is hb = 7r, (n = 1). 

The magnitude of the normal force is given by 

A(s) = 2 k 2 (s) (5 -6 k„(s)) + 3k)) m (2 R n {s) - 1). (D8) 

For the limit loop with R u m = 0, the magnitude of the normal force can be parametrized by the azimuthal angle as 
A = 2sech 4 p (5 — 6sech 2 p ). It vanishes as p —> oo, which corresponds to the points crossing the equator. 


Appendix E: Comparison of cylindrical and planar elastic curves 


To visualize the development of overhangs and the subsequent self-intersection of the loops as the length is increased, 
it is useful exploit the isometry between the cylinder and the plane. This is achieved by cutting the cylinder along 
the meridian passing through the lower turning point along the loop and “unrolling” it, as indicated in Fig. [14] If 
&m < 7t/ 2, the azimuthal extension of the loop on either side of the equator is bounded by n (A p max < 7r); see 


Figs. 14 (a)-(d). Overhangs (with p' = 0), signaling the appearance of lobes, first appear on the equator. If 
otM > 7t/2, then Ap max > 7r, and the loop necessarily develops overhangs; see Fig. |14[ e). The azimuthal extension 
grows monotonically with length: at some point Ap max = 2n and the loop makes self contact for the first time; see 
Fig. |^f). As the length is increased further, the overhangs grow indefinitely; overlapping repeatedly as they wrap 
the cylinder. 


It is also useful to compare loops on the cylinder with the corresponding planar Euler-elastic curves (dashed gray 
curves) of the same length passing through the same two points on the unrolled equator. These are superimposed on 
the unrolled loops in Fig. 14 Initially the cylindrical loop behaves like a planar elastic curve (Figs. [l4|a)-(b)). As 


the length is increased, however, the term in the potential quadratic in the normal curvature begins to play a role 
and the behavior of the two diverges, as indicated in Figs. [l4|(c)-(d). In a long loop, this term dominates: its effect is 
to elongate the cylindrical loop along the axial direction. This contrasts with its rounded planar counterpart which 
distributes its curvature equally in azimuthal and vertical directions; see[l4|e)-(f). The corresponding evolutions of the 
L 12 and L 21 are presented in Figs. 15 and 16 For the latter sequence, the cut is made at the point of self-intersection 
on the equator. In Fig. |l6jF) one can see that loops with obtuse angles extend over a azimuthal range > 27 t (the 
vertical dotted gray linesp 
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FIG. 14: Unfolding loops in the Lip sequence illustrated in Fig. [5] are represented by black solid curves, whereas their planar 
Euler-elastic counterparts (equal length and same midpoints) are represented by gray dashed curves. For small excess radius 
the two coincide, but they differ increasingly as the excess radius grows. 
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FIG. 16: Unfolded version of loops of sequence L. 2,1 illustrated in Fig. [7] are represented by black solid curves, whereas the 
planar Euler-elasticas with the same length and same midpoints are represented by gray dashed curves. 












































